Realizing acoustic qubit analogues with nonlinearly tunable phi-bits in externally driven coupled acoustic waveguides

Using experiments and theory, we investigate the behavior of nonlinear acoustic modes in a physical system composed of an array of three coupled acoustic waveguides, two of which are externally driven with two different frequencies. Nonlinear modes with frequency given by linear combinations of the driving frequencies are realizations of so-called logical phi-bits. A phi-bit is a two-state degree of freedom of an acoustic wave, which can be in a coherent superposition of states with complex amplitude coefficients, i.e., a qubit analogue. We demonstrate experimentally that phi-bit modes are supported in the array of waveguides. Using perturbation theory, we show that phi-bits may result from the intrinsic nonlinearity of the material used to couple the waveguides. We have also isolated possible effects on phi-bit states associated with the systems’ electronics, transducers and ultrasonic coupling agents used to probe the array and that may introduce extrinsic nonlinearities. These extrinsic effects are shown to be easily separable from the intrinsic behavior. The intrinsic behavior includes sharp jumps in phi-bit phases occurring over very narrow ranges of driving frequency. These jumps may also exhibit hysteretic behavior dependent on the direction of driving frequency tuning. The intrinsic states of phi-bits and multiple nonlinearly correlated phi-bits may serve as foundation for robust and practical quantum-analogue information technologies.

Quantum computing is essentially phase computing; it exploits the possibility of achieving and rotating the coherent superpositions of states of correlated multipartite systems with complex amplitudes that are represented as vectors in large, exponentially complex Hilbert spaces. The notion of "classical entanglement" for sound waves possesses the non-separability and complexity essential to reach the promise of parallelism in quantum computing, yet without the fragility of decoherence even at room temperature. In a previous publication, we have shown that one can achieve non-separable acoustic waves as classical analogues of quantum non-separable states 1 and we introduced the concept of "phi-bits" 2 . A phi-bit is a two-state degree of freedom of an acoustic wave (the acoustic spin), which can be in a coherent superposition of states with complex amplitude coefficients. Hence, a phi-bit is a qubit classical analogue-the critical component of quantum information science (QIS) platforms. Theoretically, computationally, and experimentally, we demonstrated the exponentially complex and scalable Hilbert space of states of multiple logical phi-bits (≤ 16 with 2 16 dimensional space) and the non-separability of coherent superpositions to reveal their applicability to QIS 3 . The use of phi-bits as qubit analogues for the development of quantum-analogue computing platforms that can support and navigate scalable, exponentially complex Hilbert spaces is aided by an understanding of the nonlinear phenomena that enable phi-bits as well as their correlations. Indeed, navigation of multiple correlated phi-bit Hilbert spaces to achieve operations analogous to quantum gates necessitate predictability in the phi-bit response to parametric changes in the physical system that support them. The objective of the present paper is to shed light on the origin of phi-bit behavior in a physical system composed of an externally driven array of acoustic waveguides. In particular, we focus on identifying the intrinsic nonlinearities in this physical system that lead to predictable changes in the coherent state of phi-bits upon a change in the driving conditions. We examine an array of waveguide that is driven with two different frequencies whereby phi-bits are realized as nonlinear modes with frequencies that are linear combinations of the driving frequencies. In Section "Physical system, logical phi-bits and sources of nonlinearity", we describe in detail the phi-bit-supporting physical system as well as the possible intrinsic and extrinsic www.nature.com/scientificreports/ sources of nonlinearity. Experimental results are presented in Section "Experimental results". These results focus on the behavior of the phase of the phi-bit state as one of the two driving frequency is tuned. Sections "Possible origin of the phase backgrounds" and "Possible origin of the π jumps" propose theoretical models of the phibit-supporting system that predict quantitatively and/or qualitatively experimentally observed behaviors. These models offer a window on the possible origin of sources of nonlinearity that underlie the observed behaviors.
In Section "Predictable quantum-like operations using phase jumps" we illustrate how the nonlinear response of two phi-bits may be used to operate on them. Finally in Section "Conclusions", we draw conclusions on the predictability of phi-bit behavior and their applicability to QIS.

Physical system, logical phi-bits and sources of nonlinearity
We have employed a metamaterial composed of coupled elastic wave guides to physically realize acoustic phibits 3 . The waveguides consist of three aluminum rods (6061 aluminum with diameter 1.27 cm, length 0.609 m, and density 2660 kg/m 3 ). When driven externally at two primary frequencies, an array of three elastically coupled acoustic waveguides can produce a displacement field, which when partitioned in the frequency domain leads to modes with secondary frequencies associated with logical phi-bits 3 . The waveguides are arranged in a linear array with epoxy-filled lateral gaps (50176 KwikWeld Syringe) (see Fig. 1). Ultrasonic longitudinal contact transducers (V133-RM, Olympus IMS) drive and detect the acoustic field at the rod ends. The transducers are coupled to the aluminum using commercially available honey. Separate waveform generators (B&K Precision 4055B) through PD200 amplifiers (high bandwidth, low-noise linear amplifier) excite two driving transducers at the ends of two waveguides with sinusoidal signals at frequency f 1 and f 2 . Three detecting transducers at the opposite ends collect data on the displacement field. Each detecting transducers connect to a separate channel of a Tektronix oscilloscope (MDO3024) to record signals and provide temporal and spatial information across the array.
By wrapping rubber bands (supersize bands from Walmart, 564755837) around the transducer/rod assembly and employing honey as ultrasonic coupling agent between the transducers and the rods we optimize the resolution of the acoustic modes. The rubber bands maintain a uniform pressure on the transducers throughout all experiments.
Temporal Fourier transforms of output signals generate the spectral information. The displacement field measured at the waveguide's detection-ends are the Fourier sum of modes with the primary frequencies f 1 and f 2 as well as secondary modes which frequencies are linear combination of the driving frequencies: pf 1 + qf 2 . In 3 we defined a logical phi-bit as a two-level, secondary mode of vibration whose state is characterized by the set of frequency coefficients {p, q} and spatial mode associated with two independent relative phases of the displacement between the waveguides. For an externally driven three-waveguide system, we define the state of a phi-bit, "j", by a 2 × 1 vector: 13 e iω (j) t where the angular frequency ω (j) is the linear combination of the driving angular frequencies with coefficients {p, q} . The magnitudes C 2 and C 3 are normalized to that of Picture and exploded-view of an array of three acoustic waveguides coupled with epoxy resin and schematic illustration of the experimental system for generating and detecting logical phi-bits. Three separate signal generators and amplifiers are used to drive piezoelectric transducers. Driving and detecting transducers are attached to the opposite ends of the rod-like acoustic waveguides by the pressure of three independent rubber bands. A thin layer of honey is used as ultrasonic coupling agent between the transducers and the rod ends. The signals generated by the detecting transducers enter an oscilloscope via independent input channels for analysis. The array of waveguides is suspended by thin threads for isolation.
1 are the two independent phases in waveguides 2 and 3 relative to waveguide 1. The amplitude and phases at the waveguide ends can be measured unambiguously. A single phi-bit state lives in a 2D Hilbert space h (j) and spans the Bloch sphere.
The phi-bits co-locate within the same physical space and are subject to distance-independent interactions. Tensor-product structures of systems that comprise many logical phi-bits (P) can support non-separable states in scalable, exponentially complex Hilbert spaces. Indeed, a noninteracting P phi-bit system's state is the tensor products of single phi-bit states, namely: . The tensor product of the basis vectors of single phi-bits forms a complete basis for the states of the noninteracting multi phi-bit system. This basis defines a 2 P dimensional Hilbert space H, which is the tensor product of the P Hilbert spaces of the individual noninteracting phi-bits, H = h (1) ⊗ · · · ⊗ h (P) . The secondary modes (i.e., the phi-bits) correlate via the possible interactions of the waveguide-transducer-amplifier-generator assembly. For correlated phi-bits, the Hilbert space is the same as for a noninteracting system; however, a state of the interacting system may be a separable or non-separable linear combination (with complex coefficients) of the basis vectors of H. The state of correlated phi-bits within H can be manipulated 4 to achieve quantum-like unitary operation or even algorithms that are analogous to quantum gates used in quantum computing 5,6 .
The logical phi-bits and their correlation are the result of nonlinear phenomena within the physical system of Fig. 1. Recent developments in the use of nonlinear ultrasonic techniques for non-destructive testing (NDT) of materials exploiting subharmonic and harmonic waves and guided waves 7-9 provide guidance as to what might be the factors leading to the nonlinearity in our system 10 . The support conditions of the array of waveguides may lead to nonlinearities through support/array contact. We have used suspension of the array of waveguides by two threads as it is known to minimize the possibility of support/sample nonlinearities in ultrasonic NDT 10 . It is well known that amplifiers, pre-amplifiers and signal generators exhibit so called "harmonic distortion. " That is, these electronic systems designed to generate one single desired frequency, f, also produce its harmonics 2f, 3f, etc. Indeed, we have verified that our signal generators/amplifiers exhibit some level of harmonic distortion. We therefore do not consider secondary modes that are simple harmonics of the driving frequencies f 1 or f 2 to be reliable candidate for serving as phi-bits. Secondary mode with frequency pf 1 + qf 2 may result from "intermodulation distortion". However, this commonly arises from feeding two signals with different frequencies f 1 and f 2 into a single nonlinear electronic system (signal generator and amplifiers) that will then produces the multiple frequencies af 1 + bf 2 where a, b are integers 10 . In our system, we have separated the excitations by using separate signal generator/amplifier for each applied frequency, f 1 or f 2 . In this paper, the bottom rod of the waveguide assembly is excited at f 1 and the middle rod is excited at f 2 . Separating the two excitations physically and electronically ensures that the nonlinear mixing of the frequencies occurs in the array of waveguides. However, the approach using separate amplifiers and transducers has been shown to not eliminate but only reduce the mixing of signals 10 . On the detection side, the signals generated by the three detecting transducers are acquired through three independent input channels on the oscilloscope.
Although piezoelectric materials are intrinsically nonlinear, we operate the ultrasonic transducers with amplifiers. Furthermore, we employ three individual transducers at the driving ends of the array of three waveguides to minimize the potential for frequency mixing. Fluid and gel ultrasonic coupling agents can introduce nonlinearity 10 . By applying a constant pressure at the transducer/honey/rod end interface using stretched rubber bands, we minimize the thickness of the coupling agent and reduce its potential contribution to nonlinear phenomena.
The restoring force of a stretched rubber band is intrinsically nonlinear 11 . The three rubber bands utilized to hold the driving and detecting transducers at the ends of each rod may therefore possibly contribute to the nonlinearity of the system. Finally, experimental measurements of the compression-tension asymmetry in the modulus of epoxy resins under various strain rates indicates that the epoxy used to couple the waveguides in our system may also contribute to the nonlinearities of the system 12 .

Experimental results
First, we conducted an experiment to shed light on the contribution of the rubber bands to the nonlinear behavior of the array of waveguides. We drive the system with the following frequencies: f 1 = 62 kHz and f 2 = 66 kHz . At these frequencies the wavelength of longitudinal waves in aluminum is on the order of 10 cm making propagation along the rod-like waveguides nearly one dimensional 1 . Furthermore, at these frequencies the linear longitudinal modes of the finite length waveguides are well defined and the transducers give satisfactory driving and detecting amplitudes 1 . We use a noncontact full-field scanning laser Doppler vibrometer (Polytec PSV-400) to measure the velocity field in the third waveguide as well as in the rubber band. The PSV-400 laser head was connected to a controller (Polytec OFV-5000) utilizing the VD-07 velocity decoder, which is sensitive to velocities in the range of 10 mm/s up to frequencies of MHz. Since the vibrometer can only measure velocities that are parallel to its beam, and to not rely small contraction/expansion due to Poisson's ratio perpendicular to the longitudinal direction, the laser head was placed at an oblique angle to the side of the rods to measure the longitudinal vibrations, without compromising its signal strength. Retroreflective tape was attached on the surface of the rod and rubber band to optimize the signal. The velocity time series at one point along the rod and one point along a rubber band are Fourier transformed and their corresponding frequency spectra are reported in Fig. 2. The frequency spectrum of the rod clearly shows the primary driving modes but also the presence of nonlinear modes. The modes with frequency pf 1 + qf 2 are phi-bit modes. The values of p and q are determined by the order of nonlinearity present in the system. Figure 2 shows that the current array of nonlinear waveguides enables the realization of a large number of combinations of values for p and q . The frequency spectrum of the rubber band shows little to no contribution to the nonlinear modes of vibration of the system. We can therefore neglect the contribution from www.nature.com/scientificreports/ the rubber bands to the nonlinearity of the system. All subsequent measurements of phi-bit characteristics are performed using transducers at the detection ends of the array of waveguides. Figure 3 illustrates the behavior of the experimentally measured phases ϕ 12 and ϕ 13 for the phi-bit defined by the mixed frequency 5f 1 − 2f 2 when f 1 = 62 kHz and the second driving frequency is varied in the interval [70-62 kHz] by increments of 50 Hz. This phi-bit is subsequently labelled "phi-bit A." In addition, we have measured the phases ϕ 12 f 1,2 and ϕ 13 f 1,2 for the primary modes observed at the frequencies f 1 and f 2 . Since f 1 does not change in the experiment, the associated phases remain constant. However, the frequency f 2 varying, the corresponding phases vary. We have calculated the quantities ϕ b 12 (A) = 5ϕ 12 f 1 − 2ϕ 12 f 2 and ϕ b 13 (A)= 5ϕ 13 f 1 − 2ϕ 13 f 2 which are also plotted in Fig. 3. The phases ϕ 12 (A) and ϕ 13 (A) exhibit several remarkable features in the form of upward or downward jumps superposed onto the background functions ϕ b 12 (A) and ϕ b 13 (A) . The jumps amount to phase changes on the order of π (180°). In the case of phi-bit A, there are the same number of π jumps for ϕ 12 and ϕ 13 over the range of frequency studied and they occur at the same frequencies. The jumps are very sharp and happen over ranges of frequency on the order of 100 to 200 Hz. In contrast, the calculated background phases oscillate with periods on the order of several thousands of Hz. However, other logical phi-bits exhibit different behaviors. For instance, Fig. 4 shows that the behavior of the measured phases ϕ 12 and ϕ 13 for the phi-bit defined by the mixed frequency 5f 1 − 2f 2 (labelled phi-bit "B") do not parallel each other. Here, ϕ 12 undergoes five jumps while ϕ 13 shows only one. These jumps occur over very narrow ranges of frequency not exceeding a few hundred Hz. However, similarly to phi-bit A, phi-bit B is also the superposition of jumps onto background phases ϕ b 12 (B) = 4 f 1 − 2ϕ 12 f 2 and ϕ b 13 (B)= 4ϕ 13 f 1 − 2ϕ 13 f 2 . We have also observed (not shown here) phi-bits that do not show any π jump in either ϕ 12 and ϕ 13 or phi-bits that exhibit jumps in ϕ 12 but not in ϕ 13 and vise versa. In all cases, for any phi-bit pf 1 + qf 2 , showing π jumps or not, the phases ϕ 12 and ϕ 13 superpose onto the backgrounds pϕ 12 f 1 + qϕ 12 f 2 and pϕ 13 f 1 + qϕ 13 f 2 for all frequency ranges studied. For all phi-bits, the background phases show oscillations as functions of frequency of several thousand Hz while the π jumps occur over much shorter frequency intervals of at most a few hundred Hz.
We investigated the jumps further by conducting experiments whereas the frequency f 2 is fixed at 66 kHz and f 1 is decreased by increments of 10 Hz from 59.6 kHz to 58.8 kHz or f 1 is increased by increments of 10 Hz from www.nature.com/scientificreports/ 58.8 kHz to 59.6 kHz. Figure 5 shows the difference in phase ϕ 13 for two phi-bits, namely phi-bit B and a phi-bit C with 3f 1 − 1f 2 . With a fine 10 Hz detuning, we note that the phase jumps occur over a narrow range of frequency ranging from 40 to 100 Hz. We also observe a hysteretic behavior as the frequency increases or decreases. The width of the hysteresis is less than 100 Hz. We also note that with the finer tuning interval the phase jumps are actually less than π. These jumps are closer to 160 than 180 degrees. It is worth noting that some other phi-bits we have investigated do not exhibit hysteretic phase jumps. As will be seen in Section "Possible origin of the π jumps", hysteretic behavior results from nonlinear instability which may not be present in some nonlinear modes. We have tested the system when the frequency f 1 is fixed at 62 kHz and f 2 is decreased by increments of 50 Hz from 78 kHz to 70 kHz or f 2 is increased by increments of 50 Hz from 70 kHz to 78 kHz. These experimental results (not reported here) show hysteretic phase jumps similar to those in Fig. 5. These confirm the general presence of hysteresis loops for various sets of frequencies.
These experimental results broach the challenging questions of the origin of the experimentally observed background and jumps in the phases ϕ 12 and ϕ 13 , this in the context of the sources of nonlinearity of the system. The subsequent sections develop perturbative models of nonlinear driven array of acoustic waveguides to offer possible answers to these questions.

Possible origin of the phase backgrounds
We consider a model of the acoustic metamaterial composed of three one-dimensional elastic waveguides coupled elastically along their length (Fig. 6). Each waveguide is driven externally at its end at the position x = 0 . x represents the position along the waveguides. www.nature.com/scientificreports/ The nonlinear elastic wave equation in the long wavelength limit is written as: The parameter β is proportional to the speed of sound along the waveguides. The parameter µ represents damping.
← → I is the identity matrix. α measures the elastic coupling strength between waveguides due to the epoxy. ← → M is the matrix characterizing the elastic coupling between the three waveguides. In the case of our planar array of waveguides, the coupling matrix takes the form: ⇀ F 1 and ⇀ F 2 are 3 ×1 vectors representing the external driving harmonic forces for the two different driving angular frequencies ω 1 = 2πf 1 and ω 2 = 2πf 2 .
The displacement in waveguides 1, 2 and 3 is represented by the 3 × 1 vector is a nonlinear term with strength ε . The effect of this nonlinear term will be addressed in Section "Possible origin of the π jumps". Here, we set ε =0 and solve the linear wave equation: is the linear solution. In Eq. (3), we have replaced the trigonometric functions associated with the drivers by complex exponentials. www.nature.com/scientificreports/ We can solve this equation by defining by n and ⇀ E n with n = 1, 2, 3, the eigen values and eigen vectors of the ← → M matrix, where ⇀ E n represent the spatial eigen modes across the waveguides with components E n,j , j = 1, 2, 3 . We write: The eigen modes have eigen values 1 = 0, 2 = 1 , and 3 = 3 , and are given by: We can now expand the displacement vector on the complete orthonormal basis, ⇀ E n :  Since the waveguides are finite in length, the wave numbers, k n , form a discrete set and we use a discrete summation in Eq. (9).
Equation (8) is now evaluated at x = 0 leading to the particular solutions for the driven complex amplitudes: where we define the characteristic frequency The particular linear displacement field is found to be with the complex resonant amplitudes given by Eq. (10). We note that these amplitudes are complex quantities as a result of the dissipative term iµω and therefore possess a phase. Rewriting the linear displacement field for driving forces with two different frequencies, one gets: Here, we have used two independent summation indices for the wavenumber, namely k n and k ′ n . We can rewrite the contribution of the two driving frequencies to Eq. (13) by using We can also define the displacement field at the end of the waveguides due to the coupled vibration by the renormalized 2 × 1 vectors: where C 2 and C 3 are normalized to C 1 and ϕ 12 = ϕ 2 − ϕ 1 and ϕ 13 = ϕ 3 − ϕ 1 . Here ϕ (1) 12 ,and ϕ (2) 13 correspond to the experimentally measure phases ϕ 12 f 1 , ϕ 13 f 1 , ϕ 12 f 2 , and ϕ 13 f 1 of the primary modes. These are the phases that determine the calculated phase backgrounds discussed in the experimental section. So far, our model focuses on the array of waveguides. However, in accordance with Section "Physical system, logical phibits and sources of nonlinearity", we propose a model of the array of waveguides coupled to extrinsic nonlinear damped oscillators at the rod ends that may be representative of nonlinear effects associated with the amplifiers and/or transducers and/or coupling agent: 2 (x = 0) represent the driving force on the oscillators due to the linear displacement fields of the array of waveguides at their ends, x = 0. K is some proportionality constant that converts displacement of the end of a waveguide into a driving force on the oscillator. Let us rewrite Eq. (15) for one oscillator " i": (1) i and f i and ϕ (1) i , ϕ (2) i are the i th amplitudes and phases in Eq. (14). In Eq. (16) 1 . Theses phases are linear combinations of the phase differences in the displacement fields at the ends of the waveguides in the coupled array. The linear combination is the same as that of the nonlinear frequency pω 1 + qω 2 . Since each oscillator is physically in contact with the ends of the rods, one expects that they will contribute to (15) 3 +qϕ www.nature.com/scientificreports/ the detected signals at the ends of the waveguides. The oscillators will therefore contribute a background to the phases of a phi-bit with frequency pf 1 + qf 2 that are linear combinations of the phase differences ϕ (1) 13 , ϕ (2) 12 ,and ϕ (2) 13 associated with the linear displacement of the array of waveguides. As one varies the driving frequency f 2 = ω 2 2π , the variations in ϕ (2) 12 ,and ϕ (2) 13 are controlled by the variations of the complex amplitudes A n,2 k ′ n = F (2) n ω 2 0,n k ′ n −ω 2 2 +iµω 2 as given by Eq. (10). These variations will therefore be controlled by the characteristic frequencies of the modes of vibration of the array of waveguides: ω 2 0,n k ′ n = β 2 k ′ 2 n + α 2 n . Since these characteristic frequencies are known experimentally to be spaced by on the order of several kHz between two consecutive discrete states for the same spatial mode or for states between different spatial modes 1 , we expect the background phase differences to exhibit variations with the same frequency scale. This model is therefore consistent with the observation of the background phases in the experimental measurements. It is therefore straightforward if one wants to eliminate the extrinsic effect of nonlinear oscillators coupled to the ends of the waveguides by subtracting the background phases from the overall signals.

Possible origin of the π jumps
Here, we seek analytical approximations to a variation on the nonlinear Eq. (1) with two-frequency excitation using frequency detuning with multiple time scale perturbation theory 13,14 .
Note the additional parameter δ and the dependency of damping on ε . ⇀ G( ⇀ U) models nonlinearities that act on the waveguides along their length. It may therefore represent qualitatively the effect of the epoxy nonlinearity.
We consider two time scales, τ 0 = t and τ 1 = εt . We also expend the displacement field as a sum of a zeroth order (linear) and first order (nonlinear) terms as follows: (τ 0 , τ 1 ) . The first order and second order time derivatives take the form: The wave equation to 0th order in ε is effectively the linear equation: Note that there is no damping coefficient in Eq. (25) as the effect of damping is now included into the first order equation. We will see later that the phase associated with damping included now in the 1st order equation will come back as a correction to the complex amplitude of the 0th order solution.
To 1st order, the wave equation becomes: Solutions to Eq. (25) include the sum of a solution of the homogeneous equation (i.e., without external driving) and a particular solution resulting from the driving. We have already obtained the particular solution of the 0 th order equation (linear equation) in Section "Experimental results". The complete solution takes therefore the form: For the sake of simplifying the notation, we have dropped here the different labelling of the same wave numbers k n and k ′ n used in Section "Experimental results". The first two bracket in Eq. (27a) correspond to the solution to the homogeneous equation with ω 0,n defined in Section "Experimental results". The second bracket correspond to the particular solution. Both brackets now include the complex solutions and their complex conjugate to account for the real nature of the displacement and the real form of the driving forces cosω 1 τ 0 and cosω 2 τ 0 .
We now chose a form for the term ) that enables to proceed analytically: This form assumes that the spatial modes, ⇀ E n , do not interact with each other. However, for each spatial mode, the plane wave modes, e ik n x , interact with each other. We use a third order nonlinearity for the sake of analytical tractability. The behavior of the system for even orders of nonlinearity and higher orders of nonlinearity, Q, will be discussed at the end of this section.
n=1 ⇀ E n u 1,n and using Eq. (28), we can rewrite Eq. (26) as a set of three equation, each one corresponding to a different spatial mode: To calculate u 0,n 3 , we rewrite Eq. (27b) as with η n = � n e iω 0,n τ 0 + A n,1 e iω 1 τ 0 + A n,2 e iω 2 τ 0 . We this, we rewrite the cubic term as In this expression, since we have the produce of three summations over the wave number, we reintroduced three different representations of the same wavenumbers, k n , k ′ n , k ′′ n . After several algebraic manipulation and the recognition that several multiple sums are identically the same since k n , k For this we have used the fact that dxe i(k s −k)x = δ k s ,k , the Kronecker symbol. In Eq. (33) "[c.c.]" stands for the complex conjugate of the square bracket.
Using the definition of η n and its complex conjugate, we calculate  www.nature.com/scientificreports/ These terms constitute secular terms that need to be eliminated through corresponding terms in the triple summation of Eq. (33). Seeking all the possible combinations of positive and negative k n , k ′ n , k ′′ n , that would lead to secular terms in e iω s 0,n τ 0 , we obtain the terms: We now seek the terms without e iω s 0,n τ 0 , considering for instance the terms multiplied by the exponential, e i(2ω 1 +ω 2 )τ 0 , we get the following: We introduce a detuning parameter: σ = 1 ε 2ω 1 + ω 2 − ω s 0,n , Eq. (35) is now rewritten as Adding all the secular terms in e iω s 0,n τ 0 present in the right-hand side of Eq. (33) leads to the secular equation: In Eq. (37) � ′ n = ∂� n ∂τ 1 . In the subsequent equations a prime indicates a time derivative with respect to τ 1 We now define � n (k) = 1 2 a(k)e ib(k,τ 1 ) , such that the secular equation reduces to where Ŵ 1 = 3 ω s 0,n k n A 2 n,1 + A 2 n,2 and Ŵ 2 = 1 ω s 0,n k n 2A n,1 (k s n )A n,1 (k n )A n,2 (k n ) + A n,2 (k s n )A n,1 (k n )A n,1 (k n ) . Both Ŵ 1 and Ŵ 2 are real quantities. We further define ς = σ τ 1 − b such that b ′ (k s n ,τ 1) = σ − ς ′ . We now regroup the real and imaginary terms of Eq. (38) into two separate equation and seek only steady state solution by setting all primed quantities (time derivatives with respect to τ 1 ) to zero. We finally get the set of equations: Equation (39a) and (39b) couples the amplitude of the mode of interest ( k s n ) with the amplitude of all other modes ( k n ). This equation can be simplified if we consider the self-interaction of mode k s n . In that case the last term in the square bracket reduces to 3δa 2 (k s n ) 8ω s 0,n . Note also that in considering the self interaction, one need to reduce the summation over k n in the definition of Ŵ 2 to the terms in k s n . We can get the amplitude frequency response due to self-interaction by eliminating the phase ς using the trigonometric relation: sin 2 ς + cos 2 ς = 1. That response takes the form: This result is reminiscent of the amplitude frequency response of the forced Duffing oscillator 15 . Here, we consider the features of the Duffing oscillators which are relevant to explaining the experimental behavior µiω s 0,n a k s n +iω s 0,n ib ′ k s n , τ 1 a k s n + 3 2 δa k s n k n 1 4 a 2 k s n +δa k s n Ŵ 1 ω s 0,n +δŴ 2 ω s 0,n e iσ τ 1 e −ib(k,τ 1 ) = 0 (39a) a k s n + where the amplitude decays to zero when σ = δŴ 1 . The maximum value of the amplitude is a m = |δŴ 2 | µ . We obtain the phase from Eq. (39a) as sinς = sgn(δŴ 2 ) a a m . In Fig. 7 we illustrate schematically the behavior of the amplitude and phase as functions of the frequency σ = 1 ε 2ω 1 + ω 2 − ω s 0,n . The upper overhanging side of the response leads to a jump at a m when the frequency is slowly increased. The lower overhanging side of the response is unstable, one anticipates that when the driving frequency is slowly decreased, a jump is occurring when the point J on the Fig. 7 is reached. The instability leads to a hysteresis behavior. When increasing frequency, the phase is expected to jump by π/2. Upon decreasing frequency, the phase jump will approach π. However, since the magnitude of the phase jumps will depend on the parameters δ and µ . The parameter δ which determines the strength of the nonlinearity controls the amount of overhang in the frequency domain of the amplitude. The parameter µ characterizing damping controls the width of the amplitude frequency response peak. We therefore anticipate small values of δ and µ in the case of epoxy as a source of damping and nonlinearity, leading to small overhangs and narrow amplitude and phase frequency responses.
In that case, one should observe phase jumps on the order of π with very narrow hysteresis in frequency. There may be cases of nonlinearity leading to imperceptible hysteresis.
So far, we have illustrated the effect of cubic nonlinearity on the amplitude u 0,n in forming a super-harmonic nonlinear resonance when the quantity 2ω 1 + ω 2 − ω s 0,n approaches δŴ 1 . Similar behavior may occur for other odd or even orders of nonlinearity. The mathematics, however, becomes intractable even for analytical approximate solutions. Considering an order of nonlinearity Q = p + q , in the event a nonlinear resonance occurs in the vicinity of the frequency ω s 0,n ∼ pω 1 + qω 2 , one expects phase jumps in n and therefore, u 0,n [see Eq. (27b)].
These jumps will also transfer to through Eq. (27a). Measurements of the phase differences ϕ 12 and ϕ 13 will

Predictable quantum-like operations using phase jumps
Let us now consider the application of the observed nonlinear behavior to information processing. We recall that a logical phi-bit "j" is characterized by a displacement field measured at the ends of the rods taking the general form: 13 e iω (j) t . Let us reconsider phi-bits j = A and j = B of Section "Experimental results".
The nonlinear frequencies of phi-bits A and B are 5f 1 − 2f 2 and 4f 1 − 2f 2 , respectively with f 1 = 62 kHz and f 2 = 70 kHz − �υ with �υ varying in increments of 50 Hz until it reaches 62 kHz. �υ serves therefore as a parameter to tune the state of the phi-bits. We redefine the state of the two logical phi-bits in terms of phase only by constructing the normalized state vectors: These state vectors live in single phi-bit 2D Hilbert space, with basis 1 0 , 0 1 . The quantities e iϕ (j) 12 and e iϕ (j) 13 , j = A, B , represent the complex coefficient in a linear combination of the 2D basis vectors. We now construct a state vector for a composite system composed of the two phi-bit. This state vector is represented as the tensor product: This state vector lives in the 4D Hilbert space, tensor product of two 2D single phi-bit Hilbert spaces. The basis of that 4D space is composed of the tensor products of the basis vectors of the 2D spaces. We can then define new representations of the AB phi-bit system by applying a unitary transformation. Let us consider the following block diagonal unitary matrix: Application of this transformation to � (AB) , leads to the normalized state vector  www.nature.com/scientificreports/ ← → T is one out of an infinite number of unitary transformations that produce different representations of the state of the AB composite system. Let us now see how � (AB) varies as a function of a frequency tuning parameter, �υ . In Fig. 8, we plot the phases ϕ 12 . We recall that these phases of logical phi-bits are nonlinearly correlated as they depend on the same external driving conditions. This corresponds to regions of absolute frequencies in Figs. 3 and 4 centered on 172 kHz for phi-bit A and 110 kHz for phi-bit B. We now consider a change of frequency tuning parameter from �υ 2 to �υ 1 leads to very small changes in ϕ Note that this matrix has operated on complex amplitude coefficients. By choosing that the interval [�υ 2 , �υ 1 ] encompasses an eventual hysteresis, this gate operation is reversible. This is but one example of how to realize quantum-like logic gates using logical phi-bits.
The generalization of the notion of representations of multiple (N > 2) correlated phi-bits as 2 N dimensional vector states and the exploitation of the background phases as well as nonlinear phase jumps as functions of the frequency tuning parameter may enable operations in exponentially complex Hilbert spaces.

Conclusions
We have measured experimentally and analyzed theoretically the behavior of logical phi-bits supported by a nonlinear system composed of an array of three waveguide that is driven with two different frequencies. Logical phi-bits are realized as the nonlinear modes with frequencies that are linear combinations of the driving frequencies. The coherent state of phi-bits, analogous to that of a qubit, is characterized by the phase difference between waveguides measured at the end of the waveguides. A phi-bit state lives in a 2D Hilbert space and spans the Bloch sphere. We have shown that by tuning the frequency of one of two drivers, one can change the relative phase of the waveguides and therefore navigate the phi-bit Hilbert space. The changes in phase between the waveguides exhibits a background which variations scales as thousands of Hz. This background is shown theoretically to possibly originate from extrinsic effects associated with the nonlinearity of the electronics/transducer/ultrasoniccouplant assembly that enables the driving and characterization of the physical system. This background takes the form of a linear superposition of the phases associated with the vibrational modes at the primary driving frequencies. One can therefore easily subtract this background from the experimentally measured phi-bit phases to isolate the intrinsic behavior of phi-bits. The phi-bit phases also show dramatic jumps which occur very narrow ranges of driving frequency. The phase change associated with these jumps is on the order of 180° but may vary in magnitude. These jumps also exhibit hysteresis dependent on the direction of tuning, ascending or descending, of the driving frequency. These phase variations are shown theoretically to possibly arise from intrinsic nonlinearity associated with the epoxy used to couple the waveguides along their length.
We can therefore exploit these phase jumps to achieve predictably single-bit or two-bit quantum-like operations or even quantum-like algorithms 6,16 . The phi-bits supported by the nonlinear array of waveguides are correlated. In quantum systems, correlation arises from the nonlinearity of the connection between the probability for measuring a state (the square of the wave function) and the wave function. However, since the quantum wave function is a probability amplitude, quantum computing with multiple qubits suffers from the fragility of quantum superpositions of states through the collapse of the wave function upon a perturbation such as a measurement. Here, in contrast, the intrinsic nonlinearity of the physical systems correlates the complex amplitudes (i.e., magnitude and phases) of multiple phi-bits. This classical nonlinear correlation enables us to go beyond navigating the Hilbert space of individual phi-bit by exploring multiple phi-bit superpositions of states in exponentially complex Hilbert spaces. By manipulating the state of correlated phi-bits one may be able to seek single, large unitary transformations on multiple phi-bits that can operate predictably as composites of the more conventional operations on one or two bits. The present work is therefore a step forward toward the development of decoherence-free, measurement-able, predictable, operable, correlated phi-bits-analogous to qubits-as core components for practical QIS technologies that do not suffer from quantum fragility. Future work includes establishing new correspondences between the physical state of more than two correlated logical phi-bits and representations as complex vector in high-dimensional, exponentially scaling Hilbert space. www.nature.com/scientificreports/ These representations will enable experimental implementation of a nontrivial multi phi-bit unitary operation analogous to quantum circuits by a small number of simple actions of the state of the physical system such as frequency tuning; a very challenging task for current quantum computers.

Data availability
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.